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RELATIVE MOTION OF TWO PARTICLES 
IN ELLIPTIC ORBITS 

E. R. Lancaster* 

Goddard Space Flight Center 
Greenbelt, Maryland 

NOMENCLATURE 

r = position vector at time t, r = | f | 
v = velocity vector at time t, v = | v | 

E = true anomaly at time t 
a = semimajor axis, b = l/a, c = a 1 2 
k 2 = /x = gravitational constant 
r • v = scalar product of f and v 

Over the past decade a number of papers 1 6 have presented approximate 
solutions to the problem of the relative motion of two particles in elliptic orbits 
in an inverse-square central force field. These papers have assumed that the 
relative position and velocity vectors are quite small compared to the position 
and velocity vectors of the particles, and several have further assumed one of 
the particles to be in a circular ' or nearly circular ' orbit. We derive be- 
low an exact solution to this problem, subject to no restrictions. The results 
Aerospace Engineer 
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have applications to problems of rendezvous, ncilinear error analysis, station 
keeping, targeting, surveillance, and satellite clustering. 

SOLUTION 

Kepler's equation for elliptic motion and the updating equations for position 
and velocity can be written in the form 7 

T = C + A (1 - cos C) - B sin C (1) 

r = [ 1 - H (1 - cos C)] r’ 0 + [D (1 - cos C) + N sin C] v Q (2) 

v = - (P sin C) r Q + [ 1 - S (1 - cos C)] v 0 (3) 

where 

A - r Q • v fl /k e, B = 1 - r 0 b, T = k t b/c 

H = a / r 0 , D = a (r Q • r 0 )//x, N = r„ c/k 

P = k c/ r r Q , S = a/r, C = E - E Q 

and a zero subscript indicates the value at time 0. 

Let subscript 1 on a symbol designate the value of that symbol for particle 
1 in orbit 1 and subscript 2 the value for particle 2 in orbit 2. We define 
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7 = C 2 -C 4 , e = r 2 - Tj, \ = v 2 - 


r = 


= T 2 -T,, asAj-Aj, ^ = B 2 - Bj, tj = H 2 - H, 


SrD 2 -D lf v = N 2 - Nj, p = P 2 -P,, ct = S 2 - Sj 


If we place a subscript 1 on all symbols in (1), (2), and (3) and subtract the 
resulting equations from the set with subscripts 2 on all symbols, we obtain 


T' = y + A' (1 - cos y) - B' sin y 

6 = (1 - H t F) ? 0 + (D, F + N, G) \ 0 
- (H 2 o + V F) r 20 4 (D 2 o + s F + N 2 R + v G) v 20 

\ = -P, G £ 0 + (1 - Sj F) \ 0 

- (P 2 R + p G) f 20 - (S 2 Q + a F) v 20 
where we have defined 


( 4 ) 


( 5 ) 


( 6 ) 


F = 1 - cos Cj, G = sin Cj 


(?) 


T' = r+ /9G-aF 


( 8 ) 
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O) 


A' = A 2 cos Cj + B 2 sin C, 

B' = B 2 cos Cj - A 2 sin C, 

0 = cosC 1 (l-cos>') + sinCj sin 7 

R=cosCj sin 7 - sinCj (1 - cos 7) 

To obtain equations for a, /3, r , 77 , 8 ,v,p , and a which do not suffer a loss 
of significant digits due to the subtraction of nearly equal numbers, we proceed 
as follows. 


( 10 ) 

( 11 ) 

(12) 


k (c 2 A 2 - c, Aj) = r 2Q • v 20 - r 10 • v , 0 

k [c 2 (A 2 - Aj) + A, (c 2 - c x )] = f 20 • v 20 - r 10 • v , 0 

The last equation can be solved for a. Equations for the other quantities can be 
obtained in a similar manner. The full set follows. 

k C 2 a = 7 20 '*20 “ ?io - ^,0 - k A i ( c 2 " c i> < 13 > 

ft ~ r i0 “ ^2 ^ r 20 “ r lo) (14) 
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(15) 


c 2 r = - k t (b, - b 2 ) - T, (c a - Cl ) 


r io T - a 2 - a j - H 2 (r J0 - r 10 ) 


** 8 = a 2 ( ? 20 • v 20 - 7,o • v 10 ) + (a 2 - a,) 


r io ,v io 


k ^ = C 2 ( r 20 - r io) + r lO ( C 2 - C l) 


r ,o r , P = k (c 2 - c,) - P 2 [r 2 (r 20 - r 10 ) + r , 0 (r 2 - r,)] 


r x a = a 2 - a i - S 2 («- 2 " »•,) 


f 20 ’ V 20 “ r 10 " V 10 " N) ’ r iO + C 0 ’ V 20 


r 20 “ r lO _ 6 0 ’ C r !0 + **20^ ( r iO + r 2o) 


r 2 - r , = «* (7, + 7 2 )/(r, + r 2 ) 


c 2 - c , = ( a 2 - a i) / ( C l + C 2> 


a 2 - °l = 3 1 °2 < b l - b 2> 


b l ’ b 2 = 2 ( r 20 - r ! 0 ) /r iO r 20 + V^O + V w )/m 
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In the derivation of (26) use was made of 


b = 2/r - v a /V- 


( 27 ) 


SUMMARY 


Given 7 , v 10 , 7 0 , at time 0, the steps for finding 7 and \ at time t 

follow. 

1. Compute r 20 = r J0 + e 0 , v J0 = v 10 -t k Q , r 10 = (r 1Q • r 10 ) / , v J0 = v 10 • v 10 - 

2. Compute A j . B ! » ^ * N , (equations after (3)), a j and bj from (27). 

3. Solve Kepler’s equation (1) for C, . If Cj is given rather than t, (1) is used 
to compute T t . 

4. Compute 7j from (2) and rj = (7j • 7j) 1/2 . 

5. Compute P, , Sj , Aj , B a ,H 2 , D 2 , N 2 (equations after (3)), a 2 and b 2 from 
(27). 


6. Compute in order 7 2Q -7 20 - 7, 0 *7 10 , r 20 -r 10 ,b, -b 2 ,a 2 -a,, 
c 2 - c, from (21), (22), (26), (25), (24). 

7. Compute o, t, 77, 8 , v from (13), (14), (15), (16), (17), (18). 

8. Compute F, G, T ', A', B' from (7) (8), (9), (10). 

9. Solve (4) for y . 

10. Compute Q and R from (11) and (12). 

11. Compute « from (5). 

12. Compute r* 2 = r, ■* 7 , r 2 = ( ? 2 *r 2 ) 1 2 . 
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13. Compute P 3 and S 3 (equations after (3)). 

14. Compute p and cr from (23), (19), (20). 

— • 

15. Compute \ from (6). 

Since y will usually be small, 1 - cos y should be replaced by 2 sin 2 (y/2) 
for computational purposes. 

NUMERICAL EXAMPLE 

We assume the two particles to be in coplanar circular orbits with units 
chosen such that a* = m • The initial conditions are 

x io = Vio = x io = Vio = ^ 

e„ 0 =0.001, « y0 = 0, \ k0 = 0, \ y0 = - .0004996253122 

If we let t = v /4, solve for 7j.rj.Vj . v 3 at time t and compute e and K from 
the differences r 2 - 7j and v 3 - Vj , we obtain 

e = .0015394491, 6 = -. 0001262154. 

" y 

\ = 0011853623, \ y =. 0004778069 . 

Using the formulas developed in this paper we obtain 
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€ = .001539449086, e v = - .0001262154558, 

* y 

\ = .001185362260, \ = .0004778069038. 

* y 

A ten-digit calculator was used for these calculations. 
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